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Abstract 

The inverse structure functions of exit distances have been introduced as a novel diagnostic of 
turbulence which emphasizes the more laminar regions Using Taylor's frozen field 

hypothesis, we investigate the statistical properties of the exit distances of empirical 3D fully de- 
veloped turbulence. We find that the probability density functions of exit distances at different 
velocity thresholds can be approximated by stretched exponentials with exponents varying with 
the velocity thresholds below a critical threshold. We show that the inverse structure functions 
exhibit clear extended self-similarity (ESS). The ESS exponents £(p,2) for small p (p < 3.5) are 
well captured by the prediction of £(p,2) = p/2 obtained by assuming a universal distribution of 
the exit distances, while the observed deviations for large p's characterize the dependence of these 
distributions on the velocity thresholds. By applying a box-counting multifractal analysis of the 
natural measure constructed on the time series of exit distances, we demonstrate the existence of a 
genuine multifractality, endowed in addition with negative dimensions. Performing the same anal- 
ysis of reshuffled time series with otherwise identical statistical properties for which multifractality 
is absent, we show that multifractality can be traced back to non-trivial dependence in the time 
series of exit times, suggesting a non-trivial organization of weakly-turbulent regions. 
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I. INTRODUCTION 



In isotropic turbulence, structure functions are among the favorite statistical indica- 
tors of intermittency. The (longitudinal) structure function of order p is defined by 



W S T h e K 41 theo^ obtains that W = W where e - toe average 

energy dissipation rate of the fluid element of size r and C p is a constant independent of 
Reynolds number. The K62 theory [6J] extends K41 by assuming a log-normal distribution 
of e, which was questioned by Mandelbrot |7[. The anomalous scaling properties was uncov- 
ered experimentally |J implying the non-Gaussianity of the probability distribution of the 
velocity increments. 

The velocity structure functions consider the moments of velocity increments over space. 
However, when one turns to the scalar statistics in passive scalar advection, one often con- 
siders averages of the advection time versus the distance 0, Q. An alternative quantity 



was introduced, denoted the distance structure functions 1] or inverse structure functions 

an 

T r (Sv) = (t"(Sv)), (1) 

where Sv are a set of pre-chosen thresholds of velocity increments and r(8v) is the exit 
distance defined as the minimal distance for the velocity difference to exceed Sv 

r(Sv) = inf {r : |t> i±r — Vi\ > Sv} , (2) 

given a record of velocity Vi. In the literature, alternative definitions are adopted as well, 
such as \vi +r — Vi\ > Sv or t> i+r . — Vi > Sv. 

To ensure that the exit distance is defined, the threshold Sv should be less than 5f max = 
(^max — Umin)/2, where f max and v m m are respectively the maximum and minimum of the 
record. On the other hand, there is a minimal velocity increment Sv m i n = min(|wj + i — Vi\) 
for a given record such that for any Sv < Sv m i n we have r — 1. Therefore, we consider the 
range (Sv min , Sv max ). For any Sv in this range, by construction, we will obtain finite r values 
from the velocity record. 

The statistical properties studied for synthetic data of 24630 situations from the GOY 
shell model of turbulence exhibit perfect scaling dependence of the inverse structure functions 
on the velocity threshold [jj. A completely different result was obtained in n| where an 
experimental signal was analyzed and no clear scaling was found in the exit distance structure 
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functions. For smoother stochastic fluctuations associated with a spectrum with exponent 
3 < 8 < 5, such as two-dimensional turbulence, the inverse structure functions exhibit 
bifiactafity Q. While the la rg e at fixed r of the ve.ocitv structure function emphasise 
the most intermittent region in turbulence, the large r's at fixed 5v probe the laminar regions. 
Hence, the inverse structure functions provide probes of the intermediate dissipation range 
(IDR) introduced in Q]. It is clear that the extreme events in the distribution of r 
provide the prevailing contributions to the inverse structure functions for large exponents, 
which should thus be investigated carefully. 

To our knowledge, inverse structure functions (or equivalently the statistics of exit dis- 
tances) have not been used to characterize experimental three-dimensional turbulence data. 
Here, we describe in detail the probability distribution of exit distances r and find that the 
stretched exponential distribution is a good approximation for all r's. Then, we analyze 
the convergence of the inverse structure functions and investigate their multiscaling prop- 
erties. We construct a measure based on the exit distance at each level Sv and unveil the 
multifractal nature of the measure. 



II. STANDARD PRELIMINARY TESTS ON THE EXPERIMENTAL DATA 

Very good quality high- Reynolds turbulence data have been collected at the SI ONERA 
wind tunnel by the Grenoble group from LEGI 0. We use the longitudinal velocity data 
obtained from this group. The size of the velocity time series we analyzed is N 1.73 x 10 . 

The mean velocity of the flow is approximately (v) = 20m/s (compressive effects are 
thus negligible). The root-mean-square velocity fluctuations is t> rms = 1.7m/s, leading to a 
turbulence intensity equal to / = f rms / (v) = 0.0826. This is sufficiently small to allow for 
the use of Taylor's frozen flow hypothesis. The integral scale is approximately 4m but is 
difficult to estimate precisely as the turbulent flow is neither isotropic nor homogeneous at 
these large scales. 

The Kolmogorov microscale rj is given by r\ - 
v — 1.5 x 10 _5 m 2 s _1 is the kinematic viscosity of air. dv/dt is evaluated by its discrete 
approximation with a time step increment dt = 3.5466 x 10~ 5 s corresponding to the spatial 
resolution 5 r = 0.72mm divided by (v). 

The Taylor scale is given by jl3| A = 773-7^^172 — 16.6mm. The Taylor scale is thus about 



v 2 {v) 2 
15{(dv/dt) 2 ) 



n l/4 

= 0.195mm, where 
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85 times the Kolmogorov scale. The Taylor-scale Reynolds number is Re\ = = 2000. 
This number is actually not constant along the whole data set and fluctuates by about 20%. 

We have checked that the standard scaling laws previously reported in the literature are 
recovered with this time series. In particular, we have verified the validity of the power- 
law scaling E(k) ~ with an exponent (3 very close to | over a range more than two 
decades, similar to Fig. 5.4 of [14J provided by Gagne and Marchand on a similar data set 
from the same experimental group. Similarly, we have checked carefully the determination 
of the inertial range by combining the scaling ranges of several velocity structure functions 
(see Fig. 8.6 of [iJ, Fig. 8.6]). Conservatively, we are led to a well-defined inertial range 
60 < r/rj < 2000. 

III. SCALING PROPERTIES OF INVERSE STRUCTURE FUNCTIONS 
A. The probability distributions of exit distances 

We have obtained the exit times for 26 5v values: 0.01, 0.0178, 0.0316, 0.0562, 0.1, 0.2, 0.3, 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 2.33, and 2.7144 m/s. 
Fig. ^ shows the probability density functions (pdf's) P(r) as a function of rja for different 
velocity thresholds 5v = 0.5m/s (□), 5v = lm/s (*), 5v = 1.5m/s (<) and 5v = 2m/s (*), 
where a 2 = (r 2 ) is a function of 5v. It is natural to normalize the exit distances by their 
standard deviation a(5v ) for a given Sv and obtain the pdf of these normalized exit distances 
x = r/<T as 

= aP(xa) . (3) 

The inset of Fig. plots the corresponding Q(x) for the four 5v values. One can observe 
an approximate collapse for 0.1 < x < 10 but with increasing deviations for large x's. This 
is due to the fact that the pdf's P(r) for large Sv in the semi-logarithmic plot exhibit ap- 
proximate linear behaviors over a broad range of the normalized exit distances (exponential 
distribution), while small <5v's have their pdf's with fatter tail (stretched exponential dis- 
tribution). Thus, the pdf's of exit distances are not entirely described by the single scale 
a(5v ) but are in addition slowly varying in this structure as a function of 5v . We propose 
to parameterize the shape of the pdf's as 

= Aexp[-(x/x ) m ] , for x > 1 (4) 
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FIG. 1: Empirical probability density function P(r) as a function of normalized exit distance r/a 
for 5v = 0.5, 1.0, 1.5, and 2.0 m/s. The inset shows the dependence of the corresponding <&(x) 
with respect to x = r/a defined by ©. 

where the exponent m is a function of 8v. This is quite different from the inverse statistics 
extracted from the time series of financial returns, for which the distributions of exit times 
have power law tails 3, 13 (3 Q- Actually, stretched exponential distribution is 
ubiquitous in natural and social sciences, exhibiting "fat tails" (slower decaying than expo- 
nential) with characteristic scales x [20], while power law distributions have fat tails and 
are scale-free. 

Figure El presents the fitted exponents m and the characteristic scales xq as a function of 
8v. Two regimes are observed. 

• For Sv < 1.5, the exponent m increases approximately linearly with Sv as m ~ 0.20 x 
5^ + 0.63. 

• For Sv > 1.5, m is approximately constant with a value compatible with m = 1 
corresponding to pure exponential distributions. 



B. Convergence of x p Q(x) 

A preliminary condition for analyzing the inverse structure functions is the accuracy of 
the moments of exit distances. One necessary condition for T p defined by ([TJ to converge 
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FIG. 2: The dependence of the fitted exponents m and characteristic scales xq (inset) with respect 
to 5v. 

is that the integrand r p P(r) or x p $(x) converges to zero at large r, which requires the 
closure of the integrand j2l|. We have investigated x p &(x) for different powers p and values 
5v G [0.01,2.7144] to determine how noisy is the range of r's that contribute primarily to 
T p (5v ). Fixing p, x p Q(x) is more noisy for larger Sv . For instance, x 6 $(x) clearly converges 
for 5v = 0.5m/s, but not for 5v = 2.0 m/s. We find that T p (5v) with p up to 5 can be 
evaluated with good statistical confidence. For p = 6, a reasonably good evaluation of 
T^(5v) is obtained for small 5v. The integrands for q > 8 seem divergent and the evaluation 
of the corresponding T p (5v) are less sound statistically. The typical dependence of 
as a function of x are shown in Fig. 01 for 5v = 1.5 m/s and p = 1, 2, 4, 6, 8 and 10. 

We now offer an estimate of the data size needed to estimate reliable inverse structure 
functions T p for different orders p. Let us assume that x has a stretched exponential distri- 
bution (jSJ) for x greater than some xq. Thus, the integrand of T p is I(x) = Ax p e~( x / x °^ m , 
where is a normalizing constant. We estimate that a reliable estimation of T p requires a 
good convergence of the integrant up to a value several times the value x c for which the 
integrand achieves its maximum (we use a factor k w 2 — 3 according to Fig. EJ). For the 
form of the stretched exponential distribution (0J), we have x c = xoip/m) 1 /™ 1 . On the other 
hand, the largest typical value x max observed in a sample of size iV is determined by the 
standard condition 




(5) 




FIG. 3: Plots of x p &(x) as a function of x for different values of p and velocity threshold 5v = 1.5 
m/s. 

where the integration can be performed analytically in terms of a Whittaker M function. 
When is exponential, expression (j5J) leads to the simple equation 

N = (m/A)e mx ™* . (6) 

We now write that T p (5v) is reasonably well-estimated if the range of x extends at least 
k times beyond x c . This amounts to the condition x max = kx c , where k is approximately 
independent of 5v . It follows that the minimum sample size necessary to calculate T p (Sv) 
for an exponential distribution of x's (m = 1) is given by 

N(p) = {m/A)e Kp . (7) 

For k ~ 2 as suggested from the left-middle panels of Fig. El with A ^ 1/xo with xq ~ 0.7 
(see Fig. |2I), we find N(p = 8) w 6 • 10 6 and iV(p = 9) w 5 ■ 10 7 . Thus, our data set with 
iV w 1.7 x 10 7 data points should allow us to get a reasonable estimate of the 8th order 
structure function but higher-orders become unreliable. 



C. Extended self-similarity of inverse structure functions 



To investigate the scaling properties of the inverse structure functions, we define a set of 



relative exponents using the framework of extended self-similarity (ESS) 

T p (8v) oc [T po (5v)}^\ 



221: 



(8) 



S 



where po is a reference order. In the case of velocity structure functions, po = 3 is a 
quite natural choice based on the exact Kolmogorov's four-fifth law. There is no similar 
reference for the scaling properties of the inverse structure functions and we choose somewhat 
arbitrarily p = 2. In general, ESS provides a wider scaling range for the extraction of scaling 
exponents. We will see in this subsection that Eq. (jSJ) holds for our experimental data of 
turbulence with a high accuracy. 

Figure |U presents log-log plots of T p (Sv) vs T 2 {5v) for p = 1,2, ■•• ,10 with Sv G 
[0.01,2.7144]. The straight lines hold for 0.2 < 5v < 2.7144 and over at least four or- 
ders of magnitudes in T 2 {5v ), showing the existence of extended self-similarity in the inverse 
structure functions. The scaling range for small p's seems to be broader than for large p's. 




r 2 (5v) 



FIG. 4: Dependence of the inverse structure function of order p as a function of the inverse structure 
function of order 2 taken as a reference. The straight lines exemplify the extended self-similarity 
of inverse structure functions for the threshold 5v ranging from 0.2m/s to 2.7144m/s. 

The ESS scaling exponents 2) are shown in Fig. The error bars on £(p, 2) cor- 
responds to ± one standard deviation. There is an indication that £(p, 2) has a nonlinear 
dependence as a function of p, with a downward curvature making the curve depart from 
the linear dependence 2) = p/2 observes for small p's. 

The monoscaling behavior 

£(p,2)=p/2 (9) 

is predicted from the assumption that the pdf of the normalized exit distances, and 
given by (HJ), is independent of 5v . By definition, we have T p (5v ) = (r p )g v = J °° r p Ps v (r)dr = 
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[a(Sv )] p J °° dx $$ v (x)x p . Thus, 

T p (5v) = [T 2 (5v)} 



p/2 XV 



(10) 



[/ °° dx $sv (x) x 2 ] p/2 

If <S>g v (x) is universal and independent of Sv, then the last term in (jlOj) is a number inde- 
pendent of Sv (and thus of T 2 (Sv)) and the mono-scaling £(p, 2) = p/2 follows. Thus, the 
prediction Q holds for those velocity thresholds Sv satisfying the condition that the pdf 
of exit distance is universal (independent of Sv). This is the analog of the K41 prediction 
on the standard structure functions. In our present case, there are deviations of the pdf 's 
of exit distances from the exponential law at small Sv, that we have proposed to be quan- 
tified under the form of stretched exponentials with exponents m(Sv) being a function 
of Sv as shown in Fig. |21 These deviations from exact self-similarity are weaker than for 
the direct statistics and are revealed more clearly for the higher orders. We can therefore 
attribute the deviation of the empirical 2) from the self-similarity Q at high orders to 
the non- universality of $5^(2;) which depends on the velocity levels Sv. 




FIG. 5: Dependence of the ESS exponents £(p, 2) as a function of the order p of the inverse structure 
functions. The straight line is the prediction (jHJ) obtained under the assumption of a universal pdf 
$>S v (x) of the exit distances (which is independent of the velocity thresholds 5v ). 
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IV. MULTIFRACTALITY OF THE TIME SERIES OF EXIT DISTANCES AT 
DIFFERENT 8v 



To investigate further the multifractal nature of the exit distance series {r(t) : t = 
1, ■ • • , iV = 3 x 11 x 2 19 « 1.7 x 10 7 } for a given Sv, we construct a probability measure \i 
through its integral function M(t) 

liQt 1 ,t 2 ]) = M{t 2 )-M(t 1 ) , (11) 

where M(t) = for < t < 1 and 

It] 

8=1 

for 1 < t < N. 

The box-counting method allows us to test for a possible multifractality of the measure 
jj,. The sizes s of the covering boxes are chosen such that the number of boxes of each size 
is an integer: n = N/ s G N '. We construct the partition function Z q as 

^) = X>G(*-i)^])] 9 • (13) 



i=l 



and expect it to scale as j23, 



24] 



Z q (s) ~ s <«) , (14) 



which defines the exponent r(q). For r(q), a hierarchy of generalized dimensions D q j25 



26 



27| | can be calculated according to 



D q = lim 1^ . (15) 

Do is the fractal dimension of the support of the measure. For our measure (|llj). we have 
Z^o — 1- The local singularity exponent a of the measure [i and its spectrum f(a) are related 



to r(g) through a Legendre transformation 24 1 



a(q) = dr(q)/dq , (16a) 

and 

f[a(q)] = q a (q)-r(q) . (16b) 
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We have tested the power-law scaling of Z q (s) as a function of the box size s for the 
exit time sequences at different velocity levels Sv. The scaling range is found to span over 
four orders of magnitude. Figure El plots the partition function Z q (s) for Sv = 1.5m/s as a 
function of the box size s for six different values of q in log-log coordinates. The solid lines 
are the least-square fits with power laws for each q. The correlation coefficients of the linear 
regressions (in log-log) are all larger than 0.997, demonstrating the existence of a very good 
scaling. 



FIG. 6: Plots of Zg(s) for 5v = 1.5m/s as a function of the box size s for different values of q in 
log-log coordinates. The solid lines are the least-square fits to the data using a linear regression 
(in log-log coordinates) corresponding to power laws. 

The scaling exponents r{q) are given by the slopes of the linear fits of ln[Z g (s)] as a 
function of Ins for different values of Sv. Figure [7| plots r(q) as a function of q for five 
different velocity levels Sv. The inset shows the fractal spectra f{a) obtained by the Legendre 
transformation of r(q) defined by ()16b)) . We observe that T(g)'s are concave and nonlinear, a 
diagnostic of multifractality. The maximal and minimal strength of the set of singularities, 
a max and a min , can be approximated asymptotically by lim^oo D q = lim^oo r{q)/q and 
linig^-oo D q = lim^-oo r(q)/q, respectively. It can be clearly observed that the steepness of 
the curve r(q) for negative q increases with Sv. Consequently, the maximal singularity a max 
increases with Sv, as shown in the inset of Fig. [7| where the value of a at the right endpoint 
increases with Sv. 

Since /i is conservative, r(l) = and r(0) = —1 for all Sv. For a given q < —1, the 
function r(q) decreases with Sv. For a given q > 1, there is a critical value Sv c such that 
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FIG. 7: (Color online) Scaling exponents r{q) as a function of q for different velocity levels 5v. 
The inset shows the fractal spectra f(a) obtained by the Legendre transform of r(q). 

dr(q)/d5v < when 5v < Sv c and dr(q)/dSv > when 5v > 5v c . We find that 5v c can be 
approximated by a linear function 

§ Vc = -0.058g + 0.783 , (17) 

associated with a correlation coefficient of the linear regression equal to 0.963. In addition, 
one can also see that a m i n decreases with 5v for Sv < 5v c and increases with 5v for 5v > 5v c . 
For q = —1, we find for instance that Sv c ~ 1.2m/s. 

We have performed exactly the same multifractal analysis as done above on synthetic 
time series generated from a stretched exponential distribution and on reshuffled data of the 
real exit distances. Both tests give linear scaling exponents r(q) = q — 1 in a narrower scaling 
range 64 < s < 131072, which is the earmark of monofractality. These tests strengthen the 
presence of multifractality extracted from the real exit distance data. 

In general, multifractality in time series can be attributed to either a hierarchy of changing 
broad probability density functions for the values of the time series at different scales and/or 
different long-range temporal correlations of the small and large fluctuations 2sj]. Our 
comparison, with reshuffled data and sequences with the same pdf's but no correlation 
which exhibit trivial monofractality, suggests that multifractality in the set of exit distances 
may be attributed at least in part to the existence non-trivial dependence in the time series 
of exit distances. 

An important feature of the multifractal spectrum f(a) in the inset of Fig. [7|is the ex- 

13 



321 ]. The source 



istence of negative (or latent) dimensions, that is, f(a) < [29|, |3Cj, |31 
of negative dimensions could be twofold. Firstly, the turbulent flow is a stochastic pro- 
cess, which introduces intrinsic randomness in the multifractal measure u. We note that 
negate d_ a,so appear in continue ^fractals fl Q. Se COndly , ne g a- 
tive dimensions may be interpreted geometrically by considering cuts of higher dimensional 
multifractals 

Q, H Q- 

This intuition pro pos ed by Mandelbrot has been proved math- 



ematically in the multifractal slice theorem 



371 . |38(. In the present case, the frozen 



field hypothesis is applied and we deal with one-dimensional cut of the three dimensional 
turbulence velocity field. 



V. CONCLUDING REMARKS 



Based on Taylor's frozen field hypothesis, the statistical properties of the exit distances of 
3D turbulence have been investigated. The probability density functions of exit distances at 
different velocity thresholds have been shown to be well approximated by stretched exponen- 
tials. The inverse structure functions was shown to exhibit very clear extended self-similarity 
(ESS). The ESS exponents 2) for small p < 3.5 are well described by the monofractal 
prediction 2) = p/2 obtained by assuming a universal exponential distribution of the 
exit distance. The multifractality is thus related to the dependence of the pdf's of the 
normalized exit distances on the velocity thresholds 5v. We have demonstrated that the 
sequences of exit distances for each velocity threshold Sv exhibit a clear multifractality with 
negative dimensions. The scaling ranges over which multifractality holds cover more than 
four order of magnitude in the exit distance variable. The comparison, with reshuffled data 
and sequences with the same pdf's but no correlation which exhibit trivial monofractality, 
suggests strongly that our report of multifractality is not artifactual. 

Our report of multifractality in the time series of exit distance, which tends to emphasize 
the least turbulent / most laminar regions, suggests a much richer organization of the weakly 
turbulent and close to laminar regions than believed until recently. 
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